International Journal of Engineering and Applied Sciences (IJEAS) 

ISSN: 2394-3661, Volunie-3, Issue-4, April 2016 


Missile Guidance Law Design via Backstepping 

Technique 

Hamza Zaidi, Panlong Wu, Ali Bellahcene 


Abstract — In this paper a Back-stepping Control technique 
is proposed for conunand to line-of-sight missile guidance law 
design. In this design, the three-dimensional (3-D) CLOS 
guidance problem is formulated as a tracking problem of a 
time-varying nonlinear system. Simulation results for different 
engagement scenarios illustrate the validity of the proposed 
Backstepping-based Guidance Law. 

Index Terms —Command line-of-sight (CLOS), Backstepping 
Control system, missile guidance law. 


1. Introduction 

The Concept of command to line-of-sight (CLOS) 
guidance is to oblige (force) a missile to fly as nearly as 
possible along the instantaneous line-of-sight (LOS) between 
the land tracker and the target. If the missile can continuously 
stay on the LOS, missile will intercept the target. To set 
demanded accelerations for the missile, a guidance controller 
is used at the ground station to take computation of tracker 
information about the missile and target position, angular 
velocity and acceleration of the LOS. These acceleration 
commands can then be transmitted to the missile by a radio 
link. The CLOS guidance has been identify as a low-cost 
guidance concept because it conformance placement of 
avionics on the launch platform, as opposed to mounting on 
the expendable weapons [1], [2]. Theoretically, the 

missile-target model is nonlinear and time-varying. Many 
different guidance laws have been developed over the years, 
and with the advent of highly maneuverable targets, research 
on improved guidance laws is continuing [3]-[5]. 

In this study, a Backstepping control system is proposed for 
commanding line-of-sight CLOS. The Lyapunov stability 
theorem is used to ensure the stability of the control system. 
Simulations results demonstrate the effectiveness of the 
proposed control. 

This paper is organized as follows. Formulation of 
missile-target engagement is described in Section II. The 
design procedures of the proposed Backstepping guidance 
system are constructed in Section III. Simulation results are 
set to confirm the effectiveness of the proposed control 
system in Section IV. Conclusions are drawn in Section V. 
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Table I: Definition of symbols 


Symbol 

Description 

¥, 

Yaw angle of target 

0, 

Pitch angle of target. 


Yaw angle of missile. 


Pitch angle of missile 

<l>mc 

Roll angle command. 


Azimuth angle of LOS to target. 

r, 

Elevation angle of LOS to target. 


Azimuth angle of LOS to missile 

Ym 

Elevation angle of LOS to missile 

Act 


Ay 

Ym-Yt 

% 

Gravity acceleration. 


Axial acceleration of missile 

a 

yc 

Yaw acceleration command 

a 

zc 

Pitch acceleration command. 

^yt 

Yaw acceleration of target. 

^tz 

Pitch acceleration of target. 

K 

Missile range from ground tracker. 

R, 

Target range from ground tracker. 


II. Problematique oe The Three-Dimensional CLOS 

Guidance 

The three-dimensional CLOS guidance problem shown in 
Fig. 1 is a well-known guidance model [2], which involves 
guiding the missile along the LOS to the target. The 
three-dimensional CLOS guidance model in [5, 7] will be 
repeated here for convenience. The following description in 
Table 1 will be adopted to derive the dynamic equations of 
missile. 



Fig. 1. Three-dimensional missile-target engagement diagram. 
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The origin of the inertial frame is located at the ground base. 
The axis Z/ is vertical upward and the Xj - Yj plane is 
horizontal. The origin of the missile body frame is fixed at the 
center of mass of missile, with the Xb axis forward along the 
missile centerline. The dynamics of the missile in the inertial 
frame can be represented [2] as described in (1) in the top of 
the page. 

A tracking output is defined in order to convert the CLOS 
guidance problem into a tracking problem. The LOS frame is 
shown in Fig. 2 in which the origin of the three-dimensional 
space is located at the ground base. The X^ axis forwards 
along the LOS to the missile, and the axis is horizontal to 
the left of the X^ - Tl plane. Then, the coordinates indicated in 
Fig. 2 represent the missile position in the LOS frame, and 
they are related to through rotations as follows: 


( 2 ) 


T 

The tracking output is defined as zA[zpZ 2 ] • Since and 

Z 2 cannot be measured directly, these quantities must be 

computed indirectly using the polar position data of the 
missile available from the ground tracker as 


Zi 


-^cr^ 

CG^ 

0 ~ 

_^2_ 


-SYfG^ 

-SY,SG^ 

cr,_ 


zA 


Zi 


RmdAr + r, 

_^2_ 


K AAy + y, )cy, - c(Ay + y, )sy,cAa_ 


( 3 ) 



Missile 


Ground tracker 


Ld§ _ 

Target 

Fig. 2. Definition of tracking output. 

Z 2 represents the distance from the missile to the 

LOS. Therefore, the missile will eventually intercept the 

target if the tracking output is driven to zero. The 

three-dimensional CLOS guidance problem therefore can be 
seen as a tracking problem. Define 


Note that 
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a[vi 

a[mj. 
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( 4 ) 


Using the previous equations, (1), (2), and (4) can be put 
into the following dynamic equations of missile in state-space 
form: 


x = f{x,t)+YjSj{^)‘ 

j=i 

z = h{x,t) 




where 




X, 


X, 


X, 


aM)cx,cx^ 

a,{t)sxjcx^ 




I ^..2 , ,^2 , ,, 21/2 




(x) = 


8 : 


(x) = 


(r)sx8 - g 
0 

gcxj[x;+x;+x; 

0 
0 
0 

-S<l>cmSXiCX^-C(l),^SX, 
-s^Jix^sx^ +c4>^j:x-, 

S<l>cm COSXg 

0 

0 

0 

-c(l)^j:x.,sx^+s(l)^^sx, 

-C<l>cm^X,SX^-S(l>^^CX, 

C<l>cmCXi 
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IV 4 -f V 5 -f Vg 


IV 4 + V 5 + Vg 
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^cm! (• 


2 , 2 , 2 di 

V 4 -f V 5 -f Vg 


and h{x,t) = 


hy{x,t) 


h^ix^t) 



(1) 


( 5 ) 


X^SCJ^+X^CCJ^ 


x^sy^ccj^-x^sYtSCJt+x^cYt 

The objective of CLOS guidance control is to find a control 
law to drive the tracking output z to zero. Eq. (5) can be 
rewritten as 


Zi 

Z 2 


-f 


G^^{x,t) GjjF,?) 


G2i{x,t) G22{x,t) 


u 




Fi (x, f) 

/4(x,f)_ 
zAF{x, /)+ G(x, /) • m(/) 

where Fj(x,f) = xZ^, F2{x,t)= Xl ■h2 

Gii(x,'^i * '^o *^12 -^2 * '^o *^ ’ 

G2i(x,r) = Xj • Xq-/ jj and G22(x,r)= Xj-Xg-/jj 
Xg/ii = -&,x^ccJ, - &,X2SCJ, - x^scj, + x^ca, 

= -(x,RpCr, + (y,Z2sr, - X^SCJ, + x^ca, 

XiXq/?! = -S(f)^^SX^s{X2 -CT,) + C<PcnfiX2 “ CT ,) 
^2^(h=-C<t>cmSxAXl-(xd-AmCdh -O-,) 


(5) 
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XiXoh^ = icr,cx^+sr,sx^c(x-,-( 7 ,))s<p„„+c^^^srMx 7 -cr,) 
= {cr,cx^ + sr,sx^c{Xj - a, - s(f)^^sy,s{x-i - a,) 

Xohi = ((j,sa,sr,-f,cr,sc^,)x^-{Y,cy,sa,+&,ca,sy,)x^ 
-y,x^sy, -x^ca,sy, -x^s(j,sy,+x^cy, 

= -<^,ZiSy, - y,Rp - x^ca,sy, - x^sa,sy, + x^cy, 

Xlh^ = {2a,y,sy, -a,cy, )Rp + (xfz, + (2(7,y,cy, + ^,sy,)z 2 
- 2 &,R^cy,+ 2 &,Z 2 sy, +a^{t)cx^s{x-, -a,) 

Xlhj = -{y,+&^sy,cy,)Rp-d-,ZiSy,+(&fsy,sy,+yf)z2 
- '^f,Xp - + (sxgcy, - cx^sy,c(xj - cr ,)) 

■a,{t)-gcy, 

Rp = -{<j,s<j,cy, -y,sy,ca;)Xy-{y,sy,s(T,-&,c(y,cr,)x2 
+ YXjCy, + x^cy,ca, + x^cy,sa, + x^^sy, 

= x^cy,ca, +Xscy,s<j, +x^sy, +6-,Zicy, +YZ2 
and 

Xq = 5/5f+^ y;. (x, f) • 5 / 5x,. 

i=l 
m 

7 = 1,2 

i=l 

where f-(x,t),gj - (x) and x- are the components of 
f(x,t),g(x) and x respectively. 


( 7 ) 

( 8 ) 


III. B ACKSTEPPING-BASED GUIDANCE LAW DESIGN 


Assuming that all parameters of the system ( 6 ) are known, 
the design of Backstepping control for the guidance law is 
described step-by-step as follows: 

Step 1: Define the tracking error 

(^i{t)=zAt)-z{t) ( 9 ) 

where {t) is a. desired tracking output, Then the derivative 
of tracking error can be represented as 

( 10 ) 

The z(t) can be viewed as a virtual control in above equation. 
Define the following stabilizing function 

<x(t)= i-dA)+K^e^{t) ( 11 ) 


where = 


k, 0 
0 k, 


and is a positive constant. 


The first Lyapunov function is selected as 

y,(?) = 0.5ef(r) ( 12 ) 

Step 2: Define 

e2{t) = a{t)-z(t) (13) 

Then the derivative of with respect to time is 

v; (f) = gj (fX (f) = -k^e^ (r) + gj (r)e2 (?) (14) 

Step 3: The derivative of e 2 {t )is given as 

e 2 {t) = d{t)-z{t)= z^{t)-\-k^e^{t)-F{x,t)-G{x,t)u{t) (15) 


Step 4: If all dynamics system are known, a Backstepping 
guidance law can be formulated as 

W* = G(v, t)~^ \zd (^) + (^) + (^) + ^2^2 (0 “ ^)] (1 


where K 2 


k2 0 

0 ^2 


and ^2 is a positive constant. 


Step 5: The second Lyapunov function is defined as 

V 2 (?) = Vi(?)+0.5e2(?) (17) 

Differentiating (17) and using (14) and (15), it is obtained that 

V2it) = y,it)+e2it)^2it) 

= -K^el {t)+62 {t\e^ (?) + (?) + ^jgj (?) - f(x, ?) - G{x, ?)m(?)] 

= -^,ef(?)-^ 2 e 2(?)<0 (18) 

Since > 1*^ means that gj(?) and 62 ( 1 ) 

bounded. Now define the term: 

Q(?) = K,e^ (?) + K2el (?) = - 1/2 (g, (?), g2 (?)) (19) 

then 

t 

jQ(r)?/r = y 2 (g,(o),g 2 (o))-y 2 (g,(?)ig 2 (?)) ( 20 ) 

0 

Since 1^2 (^i (o)’ ^2 (o)) is bounded and L 2 (^i(^),^ 2 (^)) i^ 
non-increasing and bounded, it can be obtained 

t 

lim [ Q(r)ir < 00 ( 21 ) 

7^00 J 

0 

AlsoQ(f)is bounded, so by using Barbalat’s Lemma [ 8 ], it 
can be shown that limQ(t) = 0 . This will imply 

?—>-CX) 

that^^(^)and^ 2 ( 0 ^^^^^^S^ as t ^ 00 .Therefore, the 

Backstepping Guidance law formulated in (16) is 
asymptotically stable. The configuration of the proposed 
Backstepping Guidance Law is shown in Fig. 3. 



Fig. 3. Backstepping Guidance System. 


IV. Numerical Simulations 

In this section, simulations are performed to illustrate the 
efficiency of the proposed Backstepping guidance law. In 
order to assess the performance characteristics in a 
closed-loop engagement scenario, it is important to specify 
target dynamics. The simplified dynamics of target motion 
can be given in the inertial frame as follows: 

x, =-ay,s\f/,-a^se,cW, 

y, =ay,cy/,-ay,s0,sii/, 

z, =a^ce,-g ( 22 ) 

W,=ay,/v,c9, 

^,={a„-gcd,)lv, 

In this paper, three simulation scenarios are examined to 
justify the effectiveness of the proposed design method. The 
simulation data and parameter data used for simulation are 
summarized in Table IT 
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Table n. Scenario and parameter data used for simulation 


States 

Scenario 1 

Scenario 2 

Scenario 3 

v(oU,(o),z,(o) [m] 

v(0U,(0),z,(0) [m/s] 

V,(0),6»,(0) [deg] 

xAo),yAo),zAo) [m] 
L(o)>L(o),4(o) [m/s] 

[deg] 

2500,5361.9,1000 

0,-340,0 

-90,0 

14.32,39.34,3.36 

70.84,151.92,28.32 

65,9.59 

5200,400,3000 

-340,0,0 

180,0 

14.56,5.43,10.01 

129.65,12.87,92.42 

20.34,32.65 

5200,5000,7500 

0,0,500 

-90,0 

28.21,34.81,26.52 

250,250,400 

45,54.73 


m/s^ 



'340 0 < f < 2 

[-44.1 t>2 


100 0<f<10 

[-44.1 f>10 


m/s^ 


/ 

< 

5g 0<t < 2.5 

-5g t > 2.5 

< 

'Og 0<f<2.5 

0.5g t > 2.5 


m/s^ 


/ 

< 

-g 0<t<2.5 

5g t > 2.5 

< 

'-g 0<r<2.5 

-g t> 2.5 

Guidance command 
frequency (Hz) 

50 

Autopilot damping 
ratio 

0.6 

Autopilot naturel 
frequency (rad/s) 

671 


<y. 


cr^ 


r. 


m 


Yt 



Fig. 4. Block diagram representation of estimation algorithm 

for guidance information. 


The first and second scenarios describes an anti-aircraft 
scenario. The third one represents an anti-missile scenario. 

Considered a 30g(g = 9.8m/^^)maneuvering limiter to limit 

the missile’s maneuverability. The pitch and yaw autopilot 
dynamics are selected to be second order linear time-invariant 
systems and the ground tracker to be a simplified differential 
tracking system with damping ratio 0.6 and nature frequency 
671 rad/s as shown in Fig. 4. The estimated values of 
cr^, and , also the measurement data of Act and Ay , 

are provided by the ground tracker. To evaluate the influence 

of measurement noise, random noises with magnitude 

2 

between ±0.3 deg are included, m/s 
The Backstepping guidance law presented in (16) is 
simulated for the same engagement scenarios. This study 
adopts the following Backstepping control law: 

= G(x, tY [za (f) + (f) +e^(t)+ ^ 2^2 (0 “ 0] 



7 0 
0 7 


and = 


20 0 
0 20 


The simulation results for scenarios 1, 2, and 3 are depicted 
in Figs.5-7, respectively. 


Table m. Miss Distance (m) 


Scenario 1 

Scenario 2 

Scenario 3 

1.8059 

2.3319 

0 .8678 


V. Conclusion 

In this paper, a Backstepping control method is applied for 
the CLOS guidance law design. Simulation results show that 
the Backstepping guidance law can achieve satisfactory 
performance and smooth missile trajectories for different 
engagement scenarios. In addition, from Table III we can 
notice those small miss distances. 
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(a) Tracking output ( Zi and Z 2 ) (a) Tracking output ( Zi and Z 2 ) 




(b) Acceleration command ( CL and ) 





(b) Acceleration command ( CL y^ and ) 



(c) Missile-target trajectory 


Fig .5 Engagement scenario 1 with Backstepping guidance law. 


Fig .6 Engagement scenario 2 with Backstepping guidance law. 
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(a) Tracking output ( Zi and Z 2 ) 



t(s) 



(b) Acceleration command (dye and d^^ ) 



Fig .7 Engagement scenario 3 with Backstepping guidance law. 
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